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ABSTRACT 

We model the star formation history (SFH) and the chemical evolution of the Galac- 
tic disk by combining an infall model and a limit-cycle model of the interstellar medium 
(ISM). Recent observations have shown that the SFH of the Galactic disk violently vari- 
ates or oscillates. We model the oscillatory SFH based on the limit-cycle behavior of the 
fractional masses of three components of the ISM. The observed period of the oscillation 
(~ 1 Gyr) is reproduced within the natural parameter range. This means that we can 
interpret the oscillatory SFH as the limit-cycle behavior of the ISM. We then test the 
chemical evolution of stars and gas in the framework of the limit-cycle model, since the 
oscillatory behavior of the SFH may cause an oscillatory evolution of the metallicity. We 
find however that the oscillatory behavior of metallicity is not prominent because the 
metallicity reflects the past integrated SFH. This indicates that the metallicity cannot 
be used to distinguish an oscillatory SFH from one without oscillations. 
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Subject headings: galaxies: ISM — Galaxy: evolution — Galaxy: stellar content - 
stars: abundances — stars: formation 

1. INTRODUCTION 

Revealing the star formation histories (SFHs) of galaxies is essential in understanding the 
galaxy formation and evolution. The SFH of the Galaxy (Milky Way) is worth studying, since 
a large number of stars are observed individually and the SFH is inferred directly from the age 
distribution of the stars. The SFH is closely related to the chemical evolution of the Galaxy. For 
example, the age-metallicity relation (e.g., Pagel 1997) of Galactic stars is generally believed to 
originate from the chemical enrichment of the Galaxy as a result of star formation. 

Eggen, Lynden-Bell, & Sandage (1962) have pioneered the modeling of the Galactic SFH and 
its chemical evolution. From the correlation between the ultraviolet excess and orbital eccentricity 
of stars, they have concluded that the Galaxy formed by collapse on a free-fall timescale from a 
single protogalactic cloud. An alternative picture of halo formation has been proposed by Searle 
& Zinn (1978). They have argued that the Galactic system formed from the capture of fragments 
such as dwarf galaxies over a longer timescale than that proposed by Eggen et al. In any case, 
determining the timescale of the infall of matter and the chemical enrichment is an important issue 
to resolve the formation history of the Galaxy. 

A number of papers have investigated the formation (e.g., Burkert, Truran, & Hensler 1992) 
and chemical evolution (e.g., Matteucci & Frangois 1989) of the Galaxy and other spiral galaxies 
(e.g., Lynden-Bell 1975; Sommer-Larsen 1996). Many models of the SFH of the Galaxy have treated 
the formation of the Galactic disk through gas infall from the halo. This scenario (the so-called 
infall model) can be consistent with the age-metallicity relation of the disk stars (e.g., Twarog 
1980), if a reasonable SFH is used. Moreover, the infall model provides a physically reasonable way 
of solving the G-dwarf problem (e.g., Pagel 1997, p. 236), contrary to the closed-box model which 
tends to overpredict the number of the low-metallicity stars. 

Since stars are formed from interstellar medium (ISM), one of the factors that determine the 
star formation rate (SFR) is the gas content of galaxies. Indeed, the SFR and the gas density 
is closely related (Kennicutt 1998). The most commonly used relation is called the Schmidt law 
(Schmidt 1959). It assumes that SFR oc p n , where p is the gas density and n = 1-2. As long as 
such a law is assumed and the infall of gas occurs continuously as a smooth function of time, the 
predicted SFH is also a smooth function of time. 

Though the "classical" (i.e., smooth) infall model is widely accepted, there are observational 
data that suggest intermittent or oscillatory star formation activities in spiral galaxies. This means 
that the SFH is not a smooth function of time. Kennicutt, Tamblyn, & Congdon (1994) have shown 
that the ratio of present-to-past SFR in spiral sample has a significant scatter. More recently, 
Tomita, Tomita & Saito (1996) have analyzed the far-infrared to .B-band flux ratio /fir//b of 
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1681 spiral galaxies (see also Devereux &; Hameed 1997). The indicator /fir//s represents the 
ratio between the present SFR and the averaged SFR over the recent Gyr. They have shown 
order-of- magnitude spread of /fir//b and suggested a violent temporal variation of the SFR. 

The intermittence of the SFH in the Galactic disk is recently suggested by Rocha-Pinto et 
al. (2000a, hereafter R00). They have provided the SFH of the Galaxy inferred from the stellar 
age of the solar neighborhood, using 552 late-type stars. The age of each star has been estimated 
from the chromospheric emission in the Can H and K lines (Soderblom, Duncan, & Johnson 1991). 
After metallicity-dependent age correction, completeness correction, and scale-height correction 2 , 
they have derived the age distribution of the stars. Then, after correcting for evolved stars, they 
have derived the SFH. They have also asserted that their SFH derived from the stars in the solar 
neighborhood is representative of the SFH in the whole disk, since the diffusion timescale of stars is 
much shorter than the Galactic age. The discussion in this paper is based on Figure 2 of R00. Based 
on their data, R00 have suggested that the star formation activity of the disk is intermittent or 
variates violently. Their suggestion has been statistically confirmed by Takeuchi & Hirashita (2000, 
hereafter TH00), who have also shown that the typical timescale of the variation is 2 Gyr. We note 
that Hernandez, Valls-Gabaud, & Gilmore (2000) have also found an oscillatory component of the 
SFH in the solar neighborhood. 

Theoretically, the intermittent or oscillatory SFH is easily reproduced if we treat the ISM as 
a nonlinear open system (Ikeuchi 1988). Ikeuchi & Tomita (1983, hereafter IT83) have considered 
the ISM composed of three phases (cold, warm, and hot) as suggested by McKee & Ostriker (1977) 
and modeled the time evolution of the fractional masses of the three components (see also Habe, 
Ikeuchi, & Tanaka 1981). Since the mass exchange among the three components is a nonlinear 
process, the limit-cycle evolution of the fractional masses can emerge (see also Scalo & Struck- 
Marcell 1986; Korchagin, Ryabstev, & Vorobyov 1994). The limit-cycle behavior is supported by 
Kamaya & Takeuchi (1997), who have interpreted the various levels of the star formation activities 
in spiral galaxies shown observationally by Tomita et al. (1996) in the framework of the nonlinear 
open system model. Their interpretation is based on the galaxy-wide limit-cycle behavior of the 
ISM. 

Another interesting topic is the chemical evolution in such a limit-cycle ISM. If the oscillatory 
SFH is considered, we may find an oscillation in a chemical enrichment process. For example, the 
age-met allicity relation of the stars in the Galactic disk may scatter because of the oscillation. We 
will examine quantitatively such a scatter caused by the limit-cycle evolution. 

In this paper, we model the oscillatory SFH in the Galactic disk proposed by R00 by combining 
the infall model and the limit-cycle model. The chemical evolution in the oscillatory SFH is also 
investigated. This paper is organized as follows. First, in § 2 we model the chemical evolution of 
the Galactic disk by using the infall model. Then, in § 3 we review the limit-cycle model of ISM. 



2 The scale height is dependent on the age of stars. 
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Some results derived from the equations are described in § 4. Finally, we discuss the SFH and the 
chemical evolution in the limit-cycle ISM in § 5. 



2. CHEMICAL EVOLUTION MODEL 

The chemical evolution model of the Galactic disk is constructed here. The model is based on 
the infall model, which is characterized by the gradual gas infall from the halo. We adopt a one- 
zone model for simplicity. In other words, the phenomena of the ISM are averaged in space. This 
simple treatment is advantageous because the response of the chemical evolution on the parameters 
is easy to examine. When we compare the result with the observational data, however, we should 
be careful whether the data are averaged or not. We comment on the one-zone approximation in 
5 3.4. 



2.1. Gas and Metal 

The changing rates of the gas mass (M g ) and metal mass (Mj, where i denotes the species of 
the metal; i = Fe, O, etc.) in the Galactic disk are described by a set of differential equations 

^ = -X,i- + E, + FX<, (2) 

where ip is the SFR, E is the total injection rate of gas from stars, F is the rate of gas infall from 
halo, Xi is the abundance of i (i.e., Xj = Mj/M g ), E± is the injection rate of element i from stars, 
and X\ is the abundance of i in the infall material (see e.g., Tinsley 1980 for the basic treatment 
of chemical evolution of galaxies). Introducing X\ enables us to treat the infall of pre-enriched gas. 
An early enrichment in the halo may be important for the initial phase of the disk formation (e.g., 
Ikuta & Arimoto 1999). 

In this paper, we choose two tracers for the metallicity, O and Fe. Almost all the oxygen is 
produced by high-mass stars, while the iron is produced mainly by Type la supernovae (SNe) as well 
as by high-mass stars. Thus, we include the contribution from Type la SNe in our formulation. We 
adopt the combination of the instantaneous recycling approximation and the delayed production 
approximation as formulated by Pagel & Tautvaisiene (1995). The evolution of Fe abundance 
based on a model of Type la SN has been considered in Kobayashi et al. (1998). With these 
approximations, E and E; L at t are expressed by using the SFR as a function of time, i(){t): 

E = R ins ^(t) + R dc Mt - t) , (3) 

Ei = [i^X i (t) + y iiiM ]^) + [12 d eiX i (t-T)+y ii dei]^(t-r), (4) 

where R and y are the returned fraction of gas from stars and the fractional mass of the newly 
formed element i, respectively, and the subscripts "ins" and "del" denote the instantaneous recycling 
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and delayed production parts, respectively. The delay time r is set to 1.3 Gyr, according to § 3 of 
Pagel & Tautvaisiene (1995). 

We assume that the disk begins to form at t = 0. Its age is assumed to be 15 Gyr. When 
t — t < 0, all the functions whose arguments depend on (t — r) are set to zero; for example, 
i/j(t — t) = when t < r. The initial condition is summarized in § 3.5. 

With the above expressions (eqs. [3] and [4]), equation (1) becomes 

^ = -(1 - R ins )m + R dc Mt -t) + F(t) , (5) 

and the combination of equations (2) and (5) leads to 
dX- 

M g ^ = Yi^W+Yi^W-^ + RwW-TftXiit-^-Xiit)] 

- Fit^X^-Xi}. (6) 



2.2. Infall Rate of Gas 

For the infall rate F, we follow TH00. They assumed an exponential form for the infall rate 

F(t) = ^exp(-t/t in ), (7) 

^in 

where Mq indicates the total mass that can fall into the galaxy. In other words, 

I Fit) dt = M . (8) 
Jo 

Normalizing equation (5) by Mq leads to 

^ = -(1 - i? ins )#) + R Ac ^{t - t) + ^- exp(-t/t in ) , (9) 

where 

™ d * s ik- <10> 

Equation (6) is also normalized by Mq as 

= Y-^^ + Y^it-^ + R^t-r^Xiit-^-X^ 

- Xl{t ]~ Xl exp(-t/t in ) . (11) 
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2.3. Star Formation Law 

In order to include the three-phase model of the ISM composed of cold, warm and hot com- 
ponents, we modify the Schmidt law with the index n = 1 (Schmidt 1959) as 

if> = M g X cold /t* , (12) 

where X co \ d is the mass ratio of the cold component to the total gas mass (M g ) and t* is the 
timescale of the cold-gas consumption to form stars. In other words, we consider that stars are 
formed from cold clouds on a gas consumption timescale of t*. The time evolution of X co i(j will be 
modeled based on IT83, which treated the ISM as a nonlinear open system, in § 3. Equation (12) 
is equivalent to 

i> = /g^coidA* • (13) 

2.4. Choice of the Parameters of the Chemical Evolution 

According to Pagel & Tautvaisiene (1995), we choose the fractional masses of newly formed 
elements (see eq. [4]) as follows: ib,ins/^Oo = 0.70, Yo,dei/^o© = 0.0, Y^W^Fe© = 0.28, and 
Ype, dei/^Fe© = 0.42, where the subscript indicates the solar value. They explained the observed 
metallicity of Galactic stars along with an infall model. For the abundances in the inflow gas, we 
examine two cases: one is the primordial case, (Xp e / Xp eQ , Xq/Xq q ) = (0, 0); the other is the 
pre-enriched case, (Xp e /X-F eQ , Xq/Xoq) = (0.1, 0.25) 3 . The parameter in the delayed production 
approximation, r, is set as r = 1.3 Gyr. The returned fractions of gas, i?i ns and Rdci, are determined 
as follows: i?j ns = 0.16, and i?dei = 0.13. The details about these values are described in Appendix 
A. We will determine t* and ii n in § 3.2. 

3. LIMIT-CYCLE MODEL OF THE ISM 

We model the oscillatory behavior of the Galactic SFH proposed by R00. IT83 have shown 
that an oscillatory behavior of the fractional masses of three components (cold, warm, and hot) 
emerges if one considers the ISM to be a nonlinear open system. An introduction and details 
concerning nonlinear open systems are found in Nicolis & Prigogine (1977). We adopt the model 
by IT83 to explain the oscillatory SFH by R00. 

As long as the infall timescale (t- m ) is much longer than the oscillatory timescale, the effect of 
the infall on the limit-cycle evolution is not significant. Indeed, as shown in Table 1, for the present 



3 This means that [Fc/H] = -1 and [O/Fe] = 0.4 in the inflow material. We note [Fe/H] = log(X Fo /Xp c0 ), 
[O/H] = log(Xo/X Q), and [O/Fe] = log(X o /Xo ) - log(X Fc /X Pc0 ). 
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case the case where i m > 9 Gyr, which is much longer than the period of oscillation proposed by 
R00 (~ 1 Gyr; see also Takeuchi & Hirashita 2000). Thus, we can apply the original model by 
IT83, which did not include the effect of infall. 



Table 1: Examined Parameters. 



Model 


t si (Gyr) 


t* (Gyr) 


tin (Gyr) 


A 


6.0 


3.9 


23 


B 


11 


7.2 


12 


C 


15 


9.8 


9 



3.1. Model Equations 

We review the model by IT83. This model is used to calculate the time evolution of the mass 
fraction of the three ISM phases (see also § 1). The result is used to calculate the SFR through 
equation (13). 

The ISM is assumed to consist of three components (McKee & Ostriker 1977); the hot rarefied 
gas (T ~ 10 6 K, n ~ 10 -3 cm -3 ), the warm gas (T ~ 10 4 K, n ~ 10 _1 cm -3 ), and cold clouds 
(T ~ 10 2 K, n ~ 10 cm -3 ). The fractional masses of the three components are A^t, A warm , and 
A co id, respectively. A trivial relation holds: 

A^hot + A warm + A cold = 1 . (14) 

The following three processes are considered (see IT83 and Ikeuchi 1988 for the details): [1] the 
sweeping of a warm gas into a cold component at the rate of aA warm (a ~ 5 x 10 -8 yr -1 ); [2] the 
evaporation of cold clouds embedded in a hot gas at the rate of 6A co idAj 2 ot (b ~ 10 -7 -10 -8 yr -1 ); 
[3] the radiative cooling of a hot gas through collisions with a warm gas at the rate of cX warm X ao t 
(c ~ 10 -6 -10~ 7 yr -1 ). Writing down the rate equations and using equation (14), we obtain 

^A co id n v v2 



dr 
dX\ lot 



= --BA co i d A hot + A(l - A co i d - A ho t) , (15) 
= -A hot (l - A cold - A ho t) + B A co id Ah ot > (16) 



dr 

where r = ct, A = a/c, and B = b/c. 

The solutions of equations (15) and (16) are classified into the following three types, according 
to the values of A and B (IT83): 

1. A > 1; all the orbits in the (Acoid, Ah G t)-plane reduce to the node (0, 1), 

2. A < 1 and B > B cr ; all the orbits reduce to a stable focus [(1 - A)/(AB + 1), A], 

3. A < 1 and B < B cr ; all the orbits converge on a limit-cycle orbit, 
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where B CT = (1 — 2A)/A 2 . Apparently, case 3 is important in the interpretation of the oscillatory 
SFH shown in R00. Thus, we choose the parameters that satisfy case 3 as will be described in the 
next subsection. 

3.2. Choice of Parameters of the Limit-Cycle Model 

Since the timescale of the variation of SFR derived by R00 is ~ 1 Gyr (see also THOO), we first 
investigate whether the period of the limit-cycle can be the order of 1 Gyr. Indeed, a Gyr-timescale 
cycle is possible in the natural parameter range. According to Figure 3 of IT83, the period can be 
~ 10 2 /c when we choose A = 0.3 and A = 0.5. Since 1/c is of the order of ~ 10 6 -10 7 yr, 10 2 /c ~ 1 
Gyr is possible. Thus, we choose A = 0.3 or A = 0.5 and 1/c = 10 7 yr to demonstrate the Gyr-scale 
oscillation of the Galactic SFH. The subsequent discussions are unchanged if we adopt another set 
of the parameters that satisfies the oscillation period of ~ 1 Gyr. 

Here we confirm that the adopted parameters are within the reasonable range of the physical 
properties of the ISM in the Galactic disk. First, Ikeuchi & Tomita (1983) estimated a from the SN 
rate and the maximum radius of an SN remnant (SNR) and obtained a ~ 5 x 10~ 8 yr -1 . Next, we 
estimate b as the reciprocal of the evaporation timescale of a cold cloud. The evaporation timescale 
estimated in Hirashita (2000a) may be applicable in the present case and we obtain b ~ 10~ 7 yr -1 . 
Finally, c is estimated from the collision rate of a cold cloud with SNRs. The collision rate t co \ is 
estimated as t co \ ~ (7ri?g NR nsNR^) _1 > where -Rsnr is the typical size of a SNR, nsNR is the number 
density of SNRs in the interstellar space, and v is the typical relative velocity between a SNR and 
a cloud. If we put i?sNR = 50 pc, nsNR = 10~ 6 pc~ 3 , 4 and v = 100 km s _1 , we obtain £ co i ~ 10 7 
yr. This means that c ~ ~ 10~ 7 yr -1 . A = 0.3 and B = 0.5 are easily satisfied if we assume 
for example a = 6 x 10 -8 yr, b = 10 -7 yr, and c = 2 x 10 -7 yr, all of which are consistent with the 
above order-of-magnitude estimates. 

THOO adopted a star formation law ip = f g /t s f and did not consider the effect of a multi- 
component medium with phase changes. In order to use their choices of the parameter values A 
and 1/c we must relate their t s { (the timescale in the classical smooth infall model) to our t* (the 
one in the oscillatory infall model). This is achieved by averaging the gas consumption rate defined 
as ip/ f g = Xcoid/t* (eq. [12]) over the whole galactic age where the oscillatory part of the efficiency 
is smoothed out and becomes l/t s f. In other words, 

$// g > = (X cold )/U = l/t rf , (17) 

where (•) indicates the time average of the quantity over the Galactic history (0 < t < 15 Gyr). 
Since X co i d = 0.65 for (A, B) = (0.3, 0.5) (Appendix B), t* = 0.65i sf . Following THOO, we examine 



4 The values of -Rsnr and wsnr are estimated according to Ikeuchi & Tomita (1983). We assume that the typical 
lifetime of a SNR is 10 7 yr. 
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the three sets of (t s f, t- m ) as summarized in Table 1. They determined the parameters by fitting 
their infall model to the observed SFH proposed by R00. All the three models provide an almost 
identical SFH and reproduce the smoothed trend of the SFH (TOO). Thus, it is meaningful to 
examine all the three cases. However, Model A gives an infall timescale much larger than that in 
e.g., Matteucci & Frangois (1989), although they assumed the same time dependence of infall as 
that in this paper. This long timescale indicates that the infall rate is almost constant over the 
history of the Galactic disk. We note that Model A predicts the highest metallicity of the three 
models (§ 4). 



3.3. Treatment of the Delayed Production 

Here we should comment on the delayed production approximation. It assumes that all 
the delayed production at t is determined by the SFR at t — r. In reality, r differs among 
Type la SNe. Thus, the delayed production is determined by the averaged SFR around t — r. 
Expecting that the lifetimes of Type la progenitors are comparable to, or longer than, 1 Gyr 
(~ the period of the SFR oscillation in the Galaxy) the averaged SFR around t — r is described as 
M g (t — r)(X co id)/t* = M g (t — r)/t s f, where M g (t — r) is the gas mass at t — r and equation (17) is 
used. Thus, we hereafter assume that 

yS(* - r) = . (18) 

tsf 



3.4. Comment on One-Zone Treatment 

Following IT83, the structure of a model galaxy is approximated by one zone. The simplicity 
of the one- zone approximation gives the advantage that the background physical processes are easy 
to see. In this subsection, we discuss the one-zone treatment. 

Habe et al. (1981) stated in their § 7 that for the one-zone assumption to be acceptable it 
is necessary that the mean distance between supernova remnants (SNRs) be less than 100 pc if a 
characteristic lifetime of SNRs of riif e ~ 10 7 yr and a mean expansion velocity of 10 km s _1 are 
adopted. This is because the SNRs should affect the whole disk for the one-zone treatment. A 
distance of less than 100 pc means that there are N ~ 10 4 SNRs in a galaxy disk, if the disk 
size of 10 kpc is assumed. This number is reasonable if SNe occur every 10 3 yr (T\\{ e /N ~ 10 7 
[yr]/10 4 ). Considering that the SN rate in a spiral galaxy is typically 1/100-1/50 yr -1 (Cappellaro 
et al. 1993), the mean distance between SNRs is less than 100 pc even if 10-20 massive stars are 
clustered in a region. 

The information of a place can travel over a distance of c Sj e fr*cyc in one period of the limit cycle, 
where c Sje s is the effective sound speed in the ISM and t cyc is the oscillatory period. Estimating 
these quantities as c Si efr ~ 10 km s" 1 (a typical velocity dispersion in the interstellar space) and 
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tcyc ~ 1 Gyr, we obtain c Sje fft cyc = 10 kpc. Thus, the information of the limit-cycle behavior can 
propagate over the whole disk. Thus, the assumption of the limit-cycle behavior over the whole 
disk may be good. In this paper, as a first step, we treat the model galaxy as being a one-zone 
object. Here we should note that the gas transport in the radial direction is difficult if we consider 
the angular momentum conservation. This difficulty may be a cause of the radial gradient of the 
metallicity and the gas-to-star fraction in spiral galaxies. 

The above discussions are not a satisfactory "proof" for the limit-cycle oscillation on the scale 
of the whole Galactic disk. (But it is important that it is not rejected.) At present, thus, the cyclic 
star formation over the whole disk is an assumption that easily explains the variation of the star 
formation activity observed in the solar neighborhood by R00. We note that it also explains the 
variety of the star formation activity of spiral galaxies (Kamaya & Takeuchi 1997). Hence, in this 
paper, we base our discussion on the limit-cycle behavior on a whole-disk scale. 

When we compare our result with the observational data, we should carefully consider to what 
extent the data is averaged. The range of Galactocentric radii that enter in the average depends 
on the age. Considering that the diffusion of stellar orbits on a scale of 1 kpc occurs in 0.2 Gyr 
(Wielen 1977), it is reasonable to assume that the observational quantities are averaged on a scale 
of more than 1 kpc. Thus, as a first step, we adopt the one-zone treatment for the Galactic disk 
to see the chemical evolution in an oscillatory SFH. We should extend our model to multi-zone 
treatment as Chiappini, Matteucci, & Gratton (1997) (see also Romano et al. 2000 for a recent 
work) in the future. Observationally, the formalism in Meusinger (1991) may be useful in order to 
link the global observed SFH with the SFH in different Galactocentric annuli. 

3.5. Initial Conditions 

The initial condition is set as follows: Xo(t = 0) = 0, XF c (t = 0) = 0, f g (t = 0) = 0, 
Acoid(* = 0) = 0.1, and Ah t(i = 0) = 0.7. The convergence to the limit-cycle occurs on the 
timescale of a few periods. The results are however not strongly dependent on the choice of the 
initial conditions for X co id and A^ot- 

4. RESULTS 

In this section, the results calculated from the equations above are presented. They are com- 
pared with the observational data. Before displaying the results, we review the solving processes 
of the equations. First, the mass fraction of the cold gas, X co \d, is calculated by equations (15) 
and (16). X co id is used to determine the SFR at t through equation (13). As for the delayed 
contribution expressed in tp(t — r), we take into account the scatter of the lifetimes of progenitors of 
Type la SNe (§ 3.3) and assume equation (18). The SFH and the chemical evolution are modeled 
by the infall model. The evolution of the gas mass normalized by the total available gas mass Mq is 
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calculated by equations (9). For t s f and t- m , we examine three cases listed in Table 1. These three 
cases are also examined in THOO. The chemical evolution is calculated by equation (11). 



4.1. Star Formation History 



The SFH calculated by our model is presented in Figure la. Since the three models predict 
almost the same SFH as indicated by THOO, we present only the result of Model A in Table 1. We 
also show the SFH observationally determined by R00 in Figure 16 in order to demonstrate the 
qualitative similarity between the model prediction and the observation. 



2.5 - 
2.0 - 

v 

^ 1.0 

0.5 
0.0 





(a) 



A 



2 4 6 8 10 12 14 
15 - t [Gyr] 



2.5 r 
2.0 : 
a 1.5: 
^ 1.0 : 

0.5 : 

0.0 ■_ 




(b) 



LrH 



2 4 6 8 10 12 14 
15 - t [Gyr] 



Fig. 1. — (a) Simulated star formation history based on our model. The star formation rate (ip) 
as a function of look-back time [(15 — t) Gyr] is normalized with the time-averaged value of ip 
((■0))- Since the three models shown in Table 1 result in almost an identical star formation history, 
only Model A is shown. (6) Star formation history derived observationally by Rocha-Pinto et al. 
(2000a). 



4.2. Metallicity Evolution 



We test the model with the metallicity data of Galactic stars. Age-metallicity relation, G- 
dwarf metallicity distribution, and [Fe/0]-[Fe/H] relation are examined. First of all, we should note 
that the yields may be uncertain because of the treatment of convection, nuclear reaction rates, 
mass loss in the asymptotic giant branch phase, etc. If the yields are systematically larger /smaller 
than assumed in this paper, the metallicities predicted by our model should be systematically 
larger /smaller. Thus, quantitative agreement by fine tuning of the parameters might be meaning- 
less. However, the qualitative behavior of the metallicity evolution in the limit-cycle ISM is not 
altered even if the yield changes. 
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4-2.1. Age-metallicity relation 

The age-metallicity relation of stars in the Galactic disk provides us with information on its 
chemical enrichment history. Thus, our model is worth testing by using the age-metallicity relation 
of the stars in the solar neighborhood. The sample is provided by Rocha-Pinto et al. (2000b), which 
used the same sample as R00. 

First, we examine the case where the infall gas is of primordial abundance [i.e., (Xp e /XF e ©, 
Xq/Xoq) = (0, 0)]. The age-metallicity relation predicted by our model is shown in Figure 2a. The 
solid, dotted, and dashed lines represent the results in Models A, B, and C, respectively. We also 
present the observational data of the age-metallicity relation by Rocha-Pinto et al. (2000b) (see their 
Table 3) . Model A predicts the highest present metallicity, since its short gas consumption timescale 
leads to the most efficient chemical enrichment. Even in Model A, however, the discrepancy between 
the model prediction and data is significant in the low-metallicity range. Rocha-Pinto et al. (2000b) 
noted the high initial metallicity of the disk and attributed it to the pre-enrichment of the gas before 
the formation of the first stars in the disk. They also have shown that the age-metallicity relation 
determined from the chromospheric age can deviate upward (i.e., metallicity is overestimated) for 
larger ages from the real relation because of the uncertainty in the age estimation. Thus, we 
also examine the case where the infalling gas in enriched with metal [(Xp e /Xp e& , Xq/Xoq) = 
(0.1, 0.25)]. The result is shown in Figure 2b. We see that the discrepancy between the model 
prediction and the observational data is reduced. From the viewpoint of modeling, the yield is also 
uncertain. Thus, we do not try any fine-tuning the age-metallicity relation. 




15 - t [Gyr] 15 - t [Gyr] 

Fig. 2. — Age-metallicity relation of the stars in the solar neighborhood. [Fe/H] is used as an 
indicator for the metallicity. The solid, dotted, and dashed lines represent the results in Models A, 
B, and C, respectively. The squares with error bars indicate the observational data point by Rocha- 
Pinto et al. (2000b). (a) Abundance in the infalling gas is assumed as (a) (A^/Ape©, Xq/Xo & ) = 
(0, 0) and (b) (XyX FcG , X^/Xoq) = (0.1, 0.25). 

In spite of such uncertainty, we can discuss the qualitative behavior of the age-metallicity 
relation. The amplitude of the oscillation of metallicity is smaller than the typical scatter of the 
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observed data points (~ 0.3 dex). The small amplitude is natural, because metallicity is determined 
by all the past history of star formation, and thus the present oscillatory star formation does not 
significantly contribute to the metallicity. 

Rocha-Pinto et al.'s data shown in Figure 2 seem to show a recent increase in [Fe/H]. The pre- 
enriched infall cannot solve this increase, since the infalling gas has too low a metallicity. However, 
we should carefully examine whether the most recent data point in Figure 2 represents the SFH of 
the whole Galactic disk, because if the recent chemical enrichment rate in the solar neighborhood is 
significantly higher than that in the whole Galaxy, the most recent data point would naturally show 
a systematically higher metallicity. Moreover, data sets shown by other authors do not necessarily 
show such an increase in the recent metallicity (e.g., Twarog 1980). It is necessary to analyze the 
observational data further before we construct a theoretical model for the recent [15 — t < 1 (Gyr)] 
increase in metallicity. 

We can expect that the oscillation behavior is more prominent for oxygen than for iron, because 
all the oxygen is produced from the "instantaneous" part (i.e., Yb,ins ^> ^b,dci)- We present the 
result for the case of pre-enriched infall [i.e., (Xp c / Xp cQ , Xq/Xo & ) = (0.1, 0.25)] in Figure 3. 
Indeed, the amplitude of the oscillation is larger than Figure 2b. However, the oscillation would 
not explain the scatter of the oxygen abundance of the stars in the Galactic disk. 




-1.5 r . . . , . . . , . . . , . . . , . . . , . . . , . . . 

2 4 6 8 10 12 14 
15 - t [Gyr] 



Fig. 3. — Age-metallicity relation of the stars in the solar neighborhood. [O/H] is used as an 
indicator for the metallicity. The solid, dotted, and dashed lines represent the results in Models A, 
B, and C, respectively. For the initial enrichment, (Xp e /Xp e& , Xq/Xq q ) = (0.1, 0.25) is assumed. 
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4-2.2. Metallicity distribution 

The G-dwarf metallicity distribution is also tested along with our model, since the primary 
motivation for the infall model is to solve the G-dwarf problem (e.g., Pagel 1997). The probability 
distribution function P(\ogX{) of the metallicity is calculated from our model as 

P(log Xi) dlog Xi = C$dt, (19) 

where the constant C is the normalization so that 

»15 Gyr 



c 



pit) (jyr 

/ 4>dt = l. (20) 

Jo 



From equation (19), we obtain the following analytical expression for P: 

P(logX t ) = C(lnlO) mx t {t) 1 , (21) 

where dXi/dt is calculated from equation (11). In comparing the distribution function with the 
observational data, we should take into account the scatter of the data. Here, we simply convolve 
P with a Gaussian kernel as 



f°° 1 

^conv(logXo) = P(u) exp 

J-oo \/2iro~ 



(log Xp - u) 
2a 2 



du , (22) 



where we adopt a = 0.1 to compare with Rocha-Pinto & Maciel (1996). We adopt these data 
because we would like to use a sample of G-dwarfs, whose lifetimes are comparable to the age of 
the universe. In Figure 4, we show P conv as a function of [Fe/H]. The solid, dotted, and dashed lines 
represent the result in Models A, B, and C, respectively. The histogram shows the data by Rocha- 
Pinto k Maciel (1996). The two figures (a) and (b) correspond to (Af c /A Fc0 , X l /X 0Q ) = (0, 0) 
and (Xp e /Xp e Q, Xq/Xoq) = (0.1, 0.25), respectively (same as Fig. 2o and b, respectively). We 
see that Model A in Figure 4b seems to be the best of all the models. However, considering the 
uncertainty in the yields, we do not try any fine tuning. The excess of the observed number of 
stars around [Fe/H] ~ 0.0 is consistent with the data by Rocha-Pinto et al. (2000b). As stated in 
§ 4.2.1, this may be due to the recent significant enrichment in the solar neighborhood. 



4.2.3. Evolution of [Fe/Oj 

In order to test the effect of the limit-cycle behavior on the [Fe/O] ratio, we examine the relation 
between [Fe/O] and [Fe/H]. Since the oxygen is mainly produced by stars with short lifetimes, the 
effect of the limit-cycle ISM is reflected by the time evolution of the oxygen abundance. On the other 
hand, the iron is also produced by stars with long lifetimes and the information of the oscillation of 
ISM phase is lost in the iron abundance. Thus, we expect that [Fe/O] oscillates as the limit-cycle 
evolution of ISM. 
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In Figure 5, we show [Fe/0]-[Fe/H] relation. The solid, dotted, and dashed lines repre- 
sent the results in Models A, B, and C, respectively. The two figures (a) and (b) correspond 
to (X f Fc /X FcQ , X { /X q) = (0, 0) and (X^/X^q, X^/Xoq) = (0.1, 0.25), respectively (same as 
Fig. 2a and b, respectively). Indeed, we see that [Fe/O] oscillates. However, the amplitude of the 
oscillation is not large. This is consistent with the age-met allicity relation (Fig. 2). As mentioned 
in § 4.2.1, the amount of metallicity reflects the past-integrated SFR; thus, as long as the mass 
of newly formed stars is not dominated in the total stellar mass, the present oscillation has little 
influence on the metallicity evolution. 

5. DISCUSSIONS 

In this paper, we have modeled the oscillatory SFH proposed observationally by R00. Our 
model is a combination of an infall model developed in the field of chemical evolution and the 
limit-cycle model proposed by Ikeuchi (1988) and his collaborators. We discuss our result in the 
following two subsections. 

5.1. Limit-Cycle Star Formation History 

The oscillatory behavior of the Galactic SFH proposed by R00 is modeled by using the limit- 
cycle model of SFH. The limit-cycle behavior of the three-phase ISM is suggested by Ikeuchi (1988) 
and his collaborators. Since the period of a limit-cycle orbit can be ~ 1 Gyr within the framework 
of Ikeuchi (1988), the Galactic SFH is explained by the limit-cycle model of the ISM. 

Recently, Hirashita & Kamaya (2000) have explained the observed scatter of star formation 
activity of a sample of spiral galaxies by using the limit-cycle model. They provided a consistent 
modeling that explains the variation in the scatter of star formation activity among the morpho- 
logical types (Sa-Sc) as shown in Kennicutt et al. (1994). Since the Galaxy is a spiral galaxy, an 
oscillatory SFH is consistent with the variation of star formation activity seen in other spirals. 

Kennicutt et al. (1994) have also presented the ratio of the present SFR to the past averaged 
SFR (indicated as b there). From Figure 2 of R00, we see that b (denoted as SFR/(SFR) in R00) 
can be as large as 2-3 for the Galactic SFH. Since the morphological type of the Galaxy is Sbc (e.g., 
Binney & Merrifield 1998, p. 171), b = 2-3 is within the range of the Sbc/Sc sample in Kennicutt 
et al. (1994, their Fig. 6). This consistency implies that the oscillatory SFH may be a common 
nature for all the spiral galaxies. 
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Fig. 4. — Distribution of the G-dwarf metallicity. The histogram shows the data by Rocha-Pinto 
k, Maciel (1996). The solid, dotted, and dashed lines represent the results in Models A, B, and C, 
respectively. All the distributions are normalized to unity when they are integrated in the whole 
range of [Fe/H]. (a) Abundance in the infalling gas is assumed as (a) (X Fe /Xp eQ , Xq/Xo q ) = (0, 0) 
and (b) (X { Fc /X Fc& , X { /X 0Q ) = (0.1, 0.25). 
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Fig. 5. — Change in [O/Fe] against [Fe/H]. The solid, dotted, and dashed lines represent the results 
in Models A, B, and C, respectively. The dot-dashed line represents the observational data as 
summarized in Fig. 3 of Kobayashi et al. (1998). (a) Abundance in the infalling gas is assumed as 
(a) (X Fc /X Fc& , X^/Xqq) = (0, 0) and (b) (X f Fc /X Fc& , X { Q /X OQ ) = (0.1, 0.25). 
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5.2. Chemical Evolution in Limit-Cycle ISM 

The chemical evolution of the Galaxy has been investigated in the framework of the limit-cycle 
ISM model. We have found that the amplitude of the oscillatory behavior of the metallicity is 
smaller than the observed scatter (Figs. 2 and 5). This indicates that the observed scatter is not 
attributed to the limit-cycle behavior. The scatter might be explained by chemical inhomogeneity 
in the Galactic disk. 

The oscillatory behavior of the metallicity is not prominent because the metallicity reflects 
all the past SFH. The integrated contribution from all the past SFH smoothes out the oscillatory 
behavior of SFR. Thus, from the viewpoint of chemical evolution we conclude that we cannot 
distinguish between the "smooth" infall model without an oscillatory behavior and the oscillatory 
infall model proposed in this paper. 

Contrary to the metallicity the dust-to-gas ratio can show a oscillatory behavior (Hirashita 
2000b). This is because the efficiency of the dust formation changes according to phase changes in 
the gas. Moreover, dust is efficiently destroyed when the mass fraction of the cold gas is small. This 
oscillation of dust amount may be important for the evolution of infrared luminosity of galaxies. 

5.3. Another Possible Mechanism for the Variation of SFR 

Rocha-Pinto et al. (2000c) gave some indication that the Magellanic Clouds could play a role 
in the SFH of the Galaxy. It is meaningful to explore their idea. Since our discussion is based on 
the limit-cycle behavior inherent in the ISM as stated in § 1, we do not include the external force 
into our formulation. Fortunately, Ikeuchi &; Tomita (1983) have investigated the behavior of the 
ISM in the presence of such an external force. Thus, we discuss the influence of the Magellanic 
Clouds based on the discussion in Ikeuchi & Tomita (1983). 

If a perturbation of the external force to the limit-cycle ISM exists, both the period and the 
amplitude of the oscillation are affected. Thus, it is not necessary that the period proper to the 
ISM is about 1 Gyr. Even in the stable-focus case (§ 3.1), an oscillation emerges. Furthermore, the 
qualitative behavior can be changed: the system can show a chaotic orbit. Since the observational 
data do not reject such a chaotic orbit affected by the perturbation of the Magellanic Clouds, the 
strongly variable SFH of the Galaxy might be due to the interaction with the Magellanic Clouds. 
However, we note that the large scatter of the star formation activity of the spiral sample (e.g., 
Tomita et al. 1996) implies a general oscillatory behavior of ISM in spiral galaxies. This is naturally 
explained if such an oscillatory behavior is caused by the limit-cycle evolution inherent in the ISM. 

Finally, we would like to note that Chiappini et al. (1997) have also considered SFR variation 
due to a different mechanism. Their rapidly variable SFR is caused by a density threshold for star 
formation: They assumed that star formation occurs only if the surface density of the gas exceeds 
a critical value, while we have assumed no threshold. However, the timescale of the SFR variation 
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in Chiappini et al. is much shorter than 1 Gyr. Although Chiappini et al.'s mechanism may indeed 
present on a timescale much shorter than 1 Gyr, it is necessary to introduce a mechanism different 
from Chiappini et al. in order to explain the SFR variability presented by R00. Thus, we have 
proposed a limit-cycle scenario for the SFR in the Galactic disk. 
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Fig. Bl. — Time evolution of the cold gas fraction Y co id (solid line). The dotted line represents 
the time-averaged value of X co id- 



A. A. DETERMINATION OF R ins AND R del 

-Ri ns and i?dci are described as 

-Rins = / (m- w m ) 4>(m) dm , (Al) 
rm, in B 

Rdel = I (jn- w m ) 4>{m) dm , (A2) 

J mi, del 

where <p(m) is the initial mass function (IMF), which is normalized so that the integral of mcf>(m) 
in the full range of the stellar mass (0.1-100M© in this paper) becomes unity; m u is the upper 
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mass cut-off of the stellar mass, and we here adopt m u = 100M©; mi iins and mi i( j e i are set as 5M 
and IMq, corresponding to the stellar lifetime of r (the parameter for the delayed production) and 
the age of galaxies. If we adopt the Salpeter's IMF (4>(m) oc m~ 2 - 35 ), we obtain R ms = 0.16 and 



The cold gas mass fraction averaged over the Galactic lifetime, (-X" C oid)i is used in § 3.2. It is 
estimated as follows. First, the time evolution of the cold gas is calculated based on the limit-cycle 
model in § 3 by adopting the parameters as A = 0.3, B = 0.5, and 1/c = 10 7 yr. The time evolution 
of X co id is shown in Figure Bl. Next, we estimate (X co \d) by averaging X co \a(t) over the Galactic 
age (Tq) as 
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